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We have investigated the polarization response of a dielectric to intense few-cycle laser pulses with 
a focus on interband tunnelling. Once charge carriers are created in an initially empty conduction 
band, they make a significant contribution to the polarization response. In particular, the coherent 
superposition of conduction- and valence-band states results in quantum beats. This quantum- 
beat part of the polarization response is affected by the excitation dynamics and attosecond-scale 
motion of charge carriers in an intense laser field. Our analysis shows that, with the onset of Bloch 
oscillations or tunnelling, the non-linear polarization response becomes sensitive to the carrier- 
envelope phase of the laser pulse. 



I. INTRODUCTION 

Multiphoton strong-field excitation of valence-band 
electrons has been intensively investigated over several 
decades. Until recently, the main motivation for exper- 
imental research has been laser damage of optical ma- 
terials (U-EI- Detailed time-resolved measurements of 
electron dynamics used to be out of reach, as the char- 
acteristic times of the most essential processes belong 
to the attosecond domain. The situation has recently 
changed due to the progress in the generation of few- 
cycle laser pulses and the development of attosecond sci- 
ence 0,H[. In this context, several recent papers should 
be pointed out. Gertsvolf et al. observed that an el- 
liptically polarized 40-fs pulse of near-infrared (NIR) ra- 
diation changes its ellipticity as it propagates through 
a few micrometres of fused silica [6j. Based on their 
numerical calculations, the authors concluded that the 
change of ellipticity is due to sub-cycle, attosecond-scale 
excitation dynamics of valence-band electrons. Their 
conclusions were confirmed by Mitrofanov et al, who 
used a non-collinear pump-probe technique to observe 
an attosecond-scale modulation of the electron density 
created by a 5-fs NIR pulse in a Si02 sample [7| . A dif- 
ferent optical strong-field phenomenon was observed by 
Ghimire et al. in ZnO exposed to high-power few-cycle 
mid-infrared laser pulses: they observed the generation of 
high harmonics with properties that cannot be explained 
by the conventional non-linear optics Q. The origin of 
these harmonics was linked to a combination of anhar- 
monic electron motion and multiple Bragg reflections at 
the boundaries of the first Brillouin zone, also known as 
"Bloch oscillations" @ . Another manifestation of Bloch 
oscillations observed by the same team was a redshift in 
the optical absorption in ZnO [lfj ■ Most recently, it was 
demonstrated that an intense few-cycle laser pulse not 
only excites charge carriers in Si02, but it also drives 
measurable currents that can be steered by controlling 
the carrier-envelope phase (CEP) [ll[ of the laser pulse 
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[12|. The theoretical interpretation of that experiment 
was given in terms of Wannier-Stark localization [l3|, [l4| 
that can be viewed as a frequency-domain description of 
Bloch oscillations. These developments created a new 
perspective for studying electron dynamics in solids ex- 
posed to strong fields, where the emphasis is placed on 
attosecond-scale dynamics. 

The foundations for the theoretical description of ex- 
tremely non-linear interaction with strong fields, where 
conventional perturbation theory breaks down, were laid 
by L. V. Keldysh in his seminal papers I15l [iff where 
multiphoton and tunnelling excitations [17] . which are 
two extreme regimes of interband transitions, were dis- 
tinguished by the Keldysh parameter 

u} Ly /mE g 
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This parameter was introduced for a monochromatic 
laser field in the approximation of parabolic bands. In 
Eq. ([1]), is the laser frequency, m = (l/m c + l/mh) -1 
is the reduced mass of an electron and a hole, E s = hujg 
is the bandgap, e is the absolute value of the electron 
charge, and Fq is the amplitude of the laser field in the 
medium. In the limit 7 » 1, where the field is weak 
and/or its frequency is large, an electron is excited from 
the valence band by absorbing an integer number of pho- 
tons in a way that is well described by conventional per- 
turbation theory. In the opposite extreme 7 < 1, where 
the field is strong and its frequency is sufficiently small, 
it is quantum tunnelling that promotes electrons to the 
conduction band. Obviously, this is the case if the exter- 
nal field is constant (ojl = 0), and interband tunnelling 
as a concept was first introduced by Zener (Til ] for time- 
independent fields. In the tunnelling regime, the exci- 
tation rate is often assumed to be determined by the 
instantaneous value of the electric field Fl (t) pj . 

Electron dynamics in the intermediate regime 7 ~ 1 
are particularly complex — a laser field cannot be treated 
as a small perturbation in this regime, and, at the same 
time, interband transitions cannot be described by pure 
tunnelling. A general expression for the cycle-averaged 
excitation rate was given by L. V. Keldysh in [l6j], but 
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a completely satisfactory description of electron dynam- 
ics within an optical cycle of an ultrashort pulse is still 
missing. The dynamics can be modelled by solving the 
time-dependent Schrodinger equation (TDSE) in a suit- 
able approximation, but there is no unambiguous way to 
distinguish between valence- and conduction-band elec- 
trons as long as the external field is present. This is 
analogous to the problem of determining the ionization 
rate of atoms and molecules in strong fields , for which 
many practical solutions exist, but there is no mathemat- 
ically rigorous solution in the case of a few-cycle pulse. 
This motivated us to study the polarization response of 
a dielectric in the regime 7 ~ 1, as interband transi- 
tions directly affect polarization and, at the same time, 
polarization is a quantum-mechanical observable that de- 
termines all optical properties of a solid. 



II. THEORY 

We are mainly interested in the strong-field regime 
where perturbation theory shows its shortcomings. In 
particular, the polarization response of a solid to a 
laser pulse is no longer described by a set of lin- 
ear and non-linear susceptibilities; therefore, quantum- 
mechanical simulations are necessary to model it. We 
solve the TDSE numerically in the velocity gauge, using 
the method described in [20| . For each value of the quasi- 
momentum k, the electron wave function is represented 
in the basis of Bloch states which are evaluated by 
solving the single-electron stationary Schrodinger equa- 
tion with an unperturbed Hamiltonian Hq: 

Ho \4>l) = E£ \K) . (2) 

Here, n is a band index. This equation is solved in the 
basis of ./V max plane waves: 

iVma,x 

<^(r)=J2^ j exp[i(k + K j )v], (3) 

where Kj denotes reciprocal lattice vectors, each of which 
can be expressed as a combination of the primitive vec- 
tors {b!,b2,b 3 }: 
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K = m a h a , m a G Z. (4) 

a=l 

Ansatz ((3]) ensures that all Bloch states satisfy <j> T ^(r + 
R) = <fik( r ) exp[ikR], where R is a vector in the direct 
Bravais lattice. This property immediately follows from 
exp[iKR] = 1. 

The coefficients in expansion ([3]) , as well as the energies 
E£, are determined by diagonalizing the Hamiltonian Hq 
in the basis of plane waves exp[iKj-r] . We assume that the 
interaction of a particular electron with lattice ions and 
other electrons is described by a pseudopotential U(r) 



and expand this pseudopotential in the plane- wave basis: 

U(v) = J2 ( 5 ) 

3=1 

U S = Tif e- iK ^C/(r)dr, (6) 

where £1 is the volume of a unit cell. Writing the unper- 
turbed Hamiltonian as 

H = -^-V 2 + U(r) (7) 
2mo 

and substituting ansatz © into Eq. @, we obtain the 
following eigensystem for the expansion coefficients C£ i 
and energies E£: 

J2 (^(k + K,) 2 % + U Ki - Kj J CJfj = £ k n q, 

(8) 

Here, mo is the electron rest mass, and the eigensystem 
can be solved independently for each quasimomentum 
k. The solutions of this eigensystem are normalized to 
satisfy 

(Cl0k) = 7^ J [CM]* 0S(r) d 3 r = S mn , (9) 

which implies Y^j |Ck j\ 2 = 1- Here and in the following, 
a scalar product (-|-) implies integration over one unit 
cell. 

Having evaluated Bloch states we use them as a 
basis to solve the TDSE and thus model the interaction 
of electrons with a laser pulse: 

|Vk(«)) - (#0 + #int W) |Vk(t)) , (10) 

where Hi nt (t) is the interaction Hamiltonian. We per- 
form our simulations in the velocity gauge, where an 
external electric field is described by the vector poten- 
tial A(r,t). We also employ the dipole approximation: 
A(r, t) = A(t). In SI units, this corresponds to using the 
following form of the interaction Hamiltonian: 

Bint (t) = — A(t)p. (11) 

m 

Note that in CGS units the right hand side of Eq. (fTl"]) 
must be divided by the vacuum speed of light c. The 
substitution of the ansatz 

JV max 

iVk(*)> = m (12) 

71=1 

into Eq. (fT0|) leads to the following system of coupled 
differential equations: 

U £^ =jES o« (t ) + J_A(t)^'rfat(*). (13) 
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Here, 



current density averaged over a unit cell: 



pi = << 



N a 



(14) 



are the matrix elements of the momentum operator p. 
An important advantage of the velocity gauge is that, as 
long as the dipole approximation holds and electron scat- 
tering is neglected, there is no coupling between different 
values of k. Mathematically, this is a consequence of 
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the Bravais lattice EL, , and we assume that both k and k' 
belong to the first Brillouin zone. Due to this property, 
Eqs. (fl"3"|) can be solved independently for each quasi- 
momentum k. According to our experience, the down- 
side of the velocity gauge is a larger number of bands 
required for convergence, as compared to length-gauge 
simulations, where different values of k are coupled to 
each other [U [22[, but even the two-band approxima- 
tion may be sufficient [23j . 

At the beginning of a simulation, all electrons are sup- 
posed to reside in the valence bands (VB), while all the 
conduction bands (CB) are empty. We solve Eqs. (p~3|) 
for each of the electrons independently and then add 
single-electron contributions together to evaluate physi- 
cal observables such as current density or polarization. 
Obviously, our approach is only applicable if (i) lat- 
tice vibrations have no significant effect on the inves- 
tigated dynamics, (ii) effects related to electron corre- 
lation and electron-hole interaction are not important, 
and (iii) distortions of the pseudopotential U(r) upon 
electronic excitations are negligible. While these approx- 
imation obviously have their limitations, single-electron 
models proved to be useful for studying basic physics 
related to phenomena occur ring on attosecond and few- 
femtosecond time scales [20I \24. |25|. 

The initial condition for solving Eqs. (TT3]) for a partic- 
ular k is 



' (t ■ ) 
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Here, n is a band where the electron was before the in- 
teraction with a laser pulse. 

In this paper, we are most interested in the macro- 
scopic electric current density j(i) induced by the laser 
pulse 



j(*)= Yl [ Jk,«Wd 3 k, 

nfVB 



(16) 



where the integral is taken over the first Brillouin zone 
an d jk. n (t) represents single-electron contributions to the 



jk,n(*) 



e 1 
m c V 



J d 3 r (Re [Vi,„(r,t)Ptfk.»(*.*)] + 
+ eA(t)\^ n (r,t)\ 2 ). (17) 



Here and in the following, we add index n to quanti- 
ties related to single-electron wave functions in order to 
indicate the initial band. Thus, jk,n(*) represents the 
contribution from all the bands to j(t) in a calculation 
where the electron initially occupied band n. The term 
containing the vector potential is required in the velocity 
gauge. 

With the aid of Eq. (|T2j) . the single-electron cur- 
rent density can be expressed via probability amplitudes 



where the summation is performed over all the vectors of a \S^)'- 



jk,„(t) = eA(t)- 



Re 



L q,l 



(<„(*) 



(18) 



Three different physical effects determine the induced 
currents: polarization due to the displacement of valence- 
band electrons, the light-driven motion of conduction- 
band electrons, and quantum beats due to coherent su- 
perpositions of valence- and conduction-band states. In 
the presence of a strong external electric field, these con- 
tributions cannot be fully separated from each other, but 
the overall induced current density j (t) is defined unam- 
biguously. It is this induced current that determines the 
change of the refractive index and the generation of new 
frequency components. However, it is conventional to de- 
scribe optical response in terms of polarization, which we 
define as 



P(i) 



j(t')dt'. 



(19) 



Note that the polarization defined this way is due to both 
bound and free electrons. Ordinarily, P(i) would be as- 
signed to bound electrons, while j(i) would describe the 
motion of free electrons. We make no attempt to make 
this distinction because, to the best of our knowledge, 
there is no rigorous way to distinguish between bound 
and free electrons as long as the external field is present 
[l9l | . It is only after the laser pulse that projecting a 
wave function onto conduction-band Bloch states yields 
physically relevant excitation probabilities. 

For simplicity, we perform our numerical simulations 
in one spatial dimension x, which we assume to be the 
polarization direction of the laser field. Most of the equa- 
tions in this section can be transformed to their one- 
dimensional forms by replacing all vector quantities and 
three-dimensional integrals with scalar quantities and 
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one-dimensional integrals: r — > x, k — ¥ k. The only ex- 
ception is Eq. (|16l) , where this procedure would result in 
j(t) being measured in wrong units. Reducing the dimen- 
sionality of our original problem, we essentially assume 
that physical observables do not depend on coordinates 
y and z, so that the integral over the first Brillouin zone 
in Eq. (|T6| reduces to 



nGVB 



jk,n(t) dk, 



(20) 



where the factor 77, measured in m~ 2 , accounts for the 
integration over k y and k z . In practice, the value of 
•q should be chosen to approximate some known opti- 
cal properties of the solid, such as its refractive index or 
absorption. 

From this point on, we will write our equations in 
atomic units (at.u.), unless stated otherwise. In atomic 
units, h — e — mo = 1, a unit of energy is equal 
to 27.21 eV, and a unit of electric field is equal to 
5.142 x 10 11 V/m. Because h = 1, we will interchangeably 
use the words "energy" and "frequency" . 

We adjusted our model potential UJx) to obtain a 
bandgap close to that of Si02: 8.95 eV [26]. To this end, 
we set the lattice constant to an = 0.5 nm = 9.45 at. u., 
which is one of the lattice constants in a-quartz, and use 
the following expression for the central unit cell: 



U(x) = -0.7 (1 + tanh[x + 0.8]) x 

x (1 + tanh[— x - 



0.8]) (21) 



for |a;| < ao/2. Outside the central unit cell, the potential 
is continued periodically: U(x + a ) = U(x). The band 
structure that corresponds to this potential is shown in 
Fig. Ha). At k — the energies E£ of the first four bands 
are -35.48, -8.33, 0.62, and 13.13 eV. We regard the first 
two bands with negative values of Eft as valence bands 
(the lowest valence band is not shown Fig. (Ha)), while 
all the bands with E% > are regarded as conduction 
bands. 

We used the following expression for the vector poten- 
tial of the (linearly polarized) laser field acting on elec- 
trons in the solid: 

A(t) = -#(t l -|;|)^cos 4 sin(w L t+^ C Ep). (22) 

Here, Fn is the amplitude of the electric field, cjl is the 
central frequency of the laser pulse, <£>cep is its carrier- 
envelope phase, 0(x) is the Heaviside step function, and 
tl is related to the full width at half maximum (FWHM) 
of \A{t)\ 2 by FWHM = 4 arccos (2" 1 / 8 ) t l /vt rs 0.5224r L . 
The FWHM of |^4(i)| 2 is very close, although not pre- 
cisely equal, to the FWHM of \F{t)\ 2 . The external elec- 
tric field acting on electrons is, by definition, 



F(t) = -A'(t). 



(23) 



Note that, in this work, we do not make any attempt to 
evaluate the screening field created by the displacement 



of electrons — this field is assumed to be a part of A(t) 
and F(t). 
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FIG. 1. (a) The uppermost valence band and the two lowest 
conduction bands of the model potential (|21|) . (b) A compar- 
ison of the refractive index evaluated according to Eq. (|24[) 
with the refractive index of fused silica [2?| . 

For a weak pulse, the polarization response is linear, 
so that, in the frequency domain, P(oj) = x( w )^( w )j 
and the refractive index can be evaluated from the linear 
susceptibility x( w ) : 



n(uj) = v/l + xH- 



(24) 



In Fig. HJb), we compare the refractive index of fused 
silica with n[uj) evaluated using our model. The re- 
fractive index is plotted against the laser wavelength 
A = 2ttc/ljJ. For these simulations, we used several pulses 
with F = 10" 5 at.u. = 5 x 10 6 V/m, FWHM = 4vr/w L , 
and values of wl that allowed us to cover the spectral 
range presented in Fig. [IJb). Normalizing the polariza- 
tion response, we set 77 = 0.111 at. u. in Eq. (|20[) and use 
this value henceforth. Given the simplicity of our model, 
we find the agreement with the measured refractive index 
very satisfactory. 

The next section reports on simulations with much 
more intense fields, where multiphoton excitations popu- 
late conduction bands. For those simulations, we had to 
use 15 bands in order to achieve numerical convergence in 
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solving the TDSE. A smaller number of bands results in 
discrepancies that are most visible in the spectral range 
occupied by low-order harmonics, while increasing the 
number of bands to 20 has a negligible effect on the po- 
larization response even for the most intense pulses that 
we used in our modelling. 



III. RESULTS AND DISCUSSION 

For all the simulations in this section, we used pulses 
with the FWHM equal to three optical oscillations: 
FWHM = 3 x 2tt/w l . The electric field F(t) of such 
a pulse is shown in Fig. [5][a) by the dotted blue line 
plotted against the number of optical cycles ui,t / (2tt) . 
The solid red line in this figure shows the polarization 
response P(t) evaluated according to Eq. (TIT)]) for the 
case where the central frequency is equal to a quarter 
of the bandgap (wl = w g /4, Al = 2ttc/lu-^ = 557 nm), 
and the peak intensity is equal to 7l = 3 x 10 13 W/cm 2 , 
which corresponds to an amplitude of the electric field of 
F = 0.029 at.u. = 1.5 V/A. The scales chosen for the 
figure emphasize that, at the leading edge of the laser 
pulse, P(t) is proportional to F(t). At later times, the po- 
larization response becomes increasingly non-linear and 
high-frequency oscillations appear, which persist even af- 
ter the laser pulse ends. The frequency of these oscilla- 
tions is close to the bandgap frequency uj g — E g /h. Thus, 
these are quantum beats appearing due to the presence 
of coherent superpositions of valence- and conduction- 
band states. This part of the polarization response will 
be the main topic of our further discussion. In order to 
distinguish it from conventional harmonics due to x^ 3 \ 
X ■> X^K ■ ■ ■ nonlinearities, we chose the central laser 
frequency to be an even divisor of the band gap, tak- 
ing advantage of the fact that harmonics of even orders 
are absent if the potential U(x) possesses the inversion 
symmetry U(— x) — U{x), like our model potential ([2"T]) . 
Note that, due to the same symmetry, the four-photon 
transition from the uppermost valence band to the low- 
est conduction band at k = is forbidden, while inter- 
band tunnelling can efficiently populate conduction-band 
states in the middle of the Brillouin zone. Indirectly, 
this is confirmed by the fact that the fast oscillations 
in Fig. Ufa) only appear at very high intensities. At a 
peak intensity of Jl = 10 13 W/cm 2 , which is just 3 times 
smaller than the one used for Fig. [21 the amplitude of 
quantum-beat oscillations is smaller by approximately a 
factor of 400, which corresponds to a drop of their inten- 
sity by five orders of magnitude. 

We investigate the quantum-beat part of the polariza- 
tion response by means of the wavelet analysis using the 
Morlet wavelet (Hi: 



W{x) 



■. exp 



2r 2 



(26) 



The parameter tw determines the size of the temporal 
window, and it must be chosen to provide an optimal 
compromise between the temporal and spectral resolu- 
tion. Our choice was tw — 3w g /cjL- 

The outcome of this analysis, S(t,uj), is shown as a 
false-colour diagram in Fig. [2jb), the bandgap energy 
being indicated by a dashed white line. For these param- 
eters, the quantum-beat signal completely dominates the 
third-order harmonic, which is barely visible in the plot. 
In this time-frequency analysis, the third-harmonic sig- 
nal is centred at 6.7 eV, and it is temporally confined to 
the FWHM of the light pulse: u} L \t\/(2n) < 1.5. 

One of the most striking features observed in Fig. [2jb) 
is the fact that the quantum-beat signal initially appears 
at frequencies exceeding cj g by roughly 2 eV/H. The in- 
stantaneous frequency of quantum beats then gradually 
decreases, approaching w g at the tailing edge of the laser 
pulse. In other words, the quantum-beat signal is neg- 
atively chirped — its lower-frequency components are de- 
layed with respect to the higher-frequency ones. This 
cannot be explained by the presence of the fifth-order 
harmonic — even though the fifth harmonic occupies the 
relevant spectral range (5/iwl = 11.2 eV), it must be con- 
fined to an even smaller temporal range than the third 
harmonic, so that S(t, u) can be considered free from 
its contribution for wit / {2-k) > 1, where most of the 
quantum-beat signal is observed. 

We seek an explanation for the chirp of the quantum- 
beat signal in the laser-driven motion of charge carriers. 
The average kinetic energy of an electron-hole pair is the 
ponderomotive energy. For a monochromatic field with 
an amplitude Fq, the ponderomotive energy is given, in 
the approximation of parabolic bands, by 



Fg 



(27) 



S(t,u) 



P(t')W(uj{t / - 1)) dt' 



(25) 



Here, m is the reduced mass of an electron and a hole. 
We obtain the effective masses of electrons and holes by 
fitting parabolas to the lowest conduction and uppermost 
valence bands in the region < n/(2ao), which yields 
m c = 0.34 at. u. and = 3.61 at. u. These values are 
in good agreement with those that can be found in lit- 
erature [23. [30I] . Thus, the reduced effective mass to in 
Eq. (J27} is equal to to = (l/m e + I/toj,)" 1 = 0.35 at.u. 
The cycle-averaged total energy of an electron-hole pair 
is equal to E g + U p (t), where U p (t) is to be evaluated by 
replacing Fq in Eq. (|27j) with the pulse envelope. This 
quantity is plotted as a red line in Fig. [2][c). It is fairly 
close to the energies where S(t, w) has local maxima with 
respect to ui. We observed this kind of negative chirp in 
numerous simulations with different parameters of the 
laser pulse, some of which are shown in Fig. [3J There- 
fore, we conclude that the frequency of quantum-beats 
exceeds w g by approximately the ponderomotive energy. 
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FIG. 2. (a) The polarization Pit) induced by a laser pulse with wl = w g /4 (Al = 557 nm), a peak intensity of 3 x 10 13 W/cm 2 , 
and ipcEP = 0. The electric field F(t) is shown by the dotted blue line; the FWHM of the pulse corresponds to the time range 
from —1.5 to 1.5 optical cycles, (b) S(t,uu), which is the result of applying the wavelet transform to P(t). The dashed white 
line marks the energy E s — hw s . (c) The solid red line is the sum of the bandgap energy Eg and the ponderomotive energy 
U p (t). For each time t, the dashed blue line depicts the above-bandgap frequency where S(t,u) has a maximum, (d) The 
quantum-beat part of the polarization response J-~ 1 [w(uj)P(uj)] (grey) and its envelope (red). The spectral window used for 
the analysis is given by Eq. (|28[) . 



Fig. E^d) presents a different kind of analysis that 
highlights the quantum-beat part of the polarization re- 
sponse. In order to suppress the third- and fifth-order 
harmonics, we multiplied the Fourier transform of P(t) 
with a soft spectral window that cuts all frequency com- 
ponents that are further than wl/2 from the bandgap 
frequency: 



w(ui) = 



2 



cos 



\oj\ — 0J i 



(28) 



V / 

The grey line in Fig.[2Jd) shows the inverse Fourier trans- 
form of w(uj)P(uj). In other words, it represents the part 
of the polarization response that oscillates at frequencies 
close to the bandgap frequency. The envelope of these 
fast oscillations is shown as a solid red line, and the dot- 
ted blue line depicts the square of the laser field. From 
this analysis, it is obvious that the quantum-beat part 
of the polarization response experiences its most rapid 
increase at a time that lies about two optical cycles af- 
ter the peak of the laser pulse. This is counter-intuitive, 



as the quantum-beat signal is due to interband excita- 
tions, and the rate of the interband excitations is ex- 
pected to reach its peak at the peak of the laser pulse. 
Another prominent feature apparent from Fig. [U[d) is 
the decay of quantum beats that begins at the trailing 
edge of the laser pulse. As our model does not account 
for any relaxation phenomena (such as electron-phonon 
interaction or spontaneous emission of radiation), the 
origin of this decay lies in the dephasing of dipole os- 
cillators associated with each pair of coherently popu- 
lated valence- and conduction-band states. Each such 
coherent superposition creates a current jk{t) described 
by Eq. (fT5|) . If the laser pulse ends at a time t\, then 



a\{t) = al(ti)exp[-i(t-h)El], so that 
j k (t) = -Re^CaJCtOr^CtOpj'e'^W-^). (29) 

q.l 

For brevity, we have omitted the index specifying the 
initial band of the electron. For a given pair of bands 
q I, Eq. (|2T)|) describes a current oscillating with a 
/c-dependent frequency E^ — E l k . Because of this depen- 
dence, the net current density (|2U|) attenuates as most 
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FIG. 3. The wavelet analysis of the polarization induced by 3-cycle pulses with a peak intensity of 3 x 10 W/cm , ipcEP = 0, and 
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currents jk(t) get out of phase with each other. 

Fig. [3] presents the wavelet analysis of polarization 
response functions for several other values of the cen- 
tral laser frequency ul- The dashed white lines depict 
E g + U p (t). In Figs.[3ja) and[3jb), these lines closely fol- 
low local maxima of S(t, uj). For smaller laser frequencies, 
the ponderomotive potential U p signihcantly exceeds the 
width of the lowest conduction band (see Fig.QJa)). This 
is the regime where Bloch oscillations [3l[ must play an 
important role Q. 

In Fig. [3Ja) , where valence-band electrons can be ex- 
cited by five-photon absorption, we made an exception 
from our choice to use those values of Wl that are even 
divisors of uj g . Still, the time- frequency analysis looks 
very similar to that for ojl = w g/4, shown in Fig. [Hh), 
and the contribution from the fifth harmonic generation 
is so weak that it is practically invisible in our colour 
scheme. For all the cases shown in Fig. |3l we observe 
the same qualitative features as those in Fig. &b): the 
frequency components of the polarization response that 
lie above the bandgap appear delayed with respect to the 
peak of the laser pulse, and the instantaneous frequency 
of these oscillations decreases with time, eventually ap- 
proaching Lu g . We also observe a rapid decrease in the 
magnitude of S(t,Lu). 

To further understand the origin of these trends, we 



plot the sum of all conduction-band populations at the 
end of our simulations in Fig. @] These are the probabil- 
ities to find an electron with a given quasimomentum k 
in one of the conduction bands at a time t > tl where 
the laser field is absent (since we do not account for re- 
laxation phenomena, populations do not change in the 
absence of external fields). The distribution of excited 
electrons dramatically changes as ojl decreases from uj g /4 
to ojg/10: the pronounced peak at k = disappears, and 
many irregular peaks appear across the Brillouin zone. 
At the same time, the distributions become increasingly 
sensitive to the carrier-envelope phase of the laser pulse. 
We consider two physical phenomena that may be re- 
sponsible for this trend: interband tunnelling and Bloch 
oscillations. 

The first phenomenon, interband tunnelling, repre- 
sents an extreme regime of interband excitations where 
the Keldysh parameter 7 is much smaller than 1. In this 
regime, interband transitions should be viewed not as a 
result of absorbing a certain number of photons, but as a 
result of quantum tunnelling [18j. Also, in the tunnelling 
regime, the carrier-envelope phase is an important pa- 
rameter because it is the field of the laser pulse, rather 
than its envelope, that controls the tunnelling rate. Ac- 
cording to Fig. 3] and all the other similar simulations 
that we have performed, the onset of the CEP depen- 
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dence coincides with 7 decreasing below a value approx- 
imately equal to 0.7. Furthermore, the tunnelling pic- 
ture provides an intuitive explanation for the multitude 
of peaks in Figs.Hfc) and0Jd): they may be attributed to 
the interference among electron wave packets launched to 
conduction bands within different half-cycles of the laser 
pulse, which is analogous to the above-threshold ioniza- 
tion spectra of atoms exposed to intense laser pulses [HJ . 



What the tunnelling picture alone cannot explain is 
the disappearance of the peak at k = 0. Indeed, all 
analytical models for the tunnelling rate predict that it 
should rapidly decrease as the bandgap increases 0, [H, 
[33|], and our model potential has the smallest bandgap 
at k — 0. Therefore, the motion of excited electrons in 
conduction bands must play a significant role 23] , Bloch 
oscillations being the most prominent manifestation of 
this motion. An electron excited to the bottom of a band 
at a zero-crossing of the vector potential (A(to) = 0) 
may reach the boundary of the first Brillouin zone if, at 
some later moment t\, the condition 1^4(^)1 = n/a is 
satisfied. Thus, Bloch oscillations should be considered 
important if the amplitude of the vector potential Fq /wl 



is comparable to 7r/a . Introducing the Bloch frequency 
w B = (SI units), (30) 

this condition can be expressed as ub/wl = tt- The la- 
bels in Figs. SJa-d) provide values of ubM in units of 
7r for all the four simulations. The appearance of the 
CEP-dependence and the complex structures in the dis- 
tribution of conduction-band electrons coincide not only 
with 7 becoming smaller than 1, but also with ub/wl be- 
coming larger than 7r. In fact, this is not a coincidence, 
as 

^ = (SI units), (31) 

as follows from Eqs. ([TJ and (J30j) . For our parameters, 
wb/wl ~ 1.027r/7, so that the onset of tunnelling coin- 
cides with the onset of Bloch oscillations, irrespectively 
of the laser pulse parameters. 

Without intcrband excitations, Bloch oscillations ob- 
viously cannot occur, but is it the nature of interband 
transitions or the Bloch oscillations themselves that are 
mainly responsible for the observed CEP dependencies? 
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While we cannot answer this question with certainty, we 
can argue that interband tunnelling alone is not sufficient 
for a polarization response to be CEP-dependent. In- 
deed, tunnelling must be the dominant excitation mech- 
anism in the parameter regime of Fig. |4ja) because, in 
spite of the forbidden four-photon transitions at k = 0, 
the bottom of the conduction band is strongly populated 
after the laser pulse. Yet the polarization response is al- 
most CEP-independent. Therefore, our simulations with 
wl = provide an example where the dominance of 

interband tunnelling over (perturbative) multiphoton ex- 
citations does not lead to a significant CEP dependence. 

The quantum-mechanical observables that we have dis- 
cussed so far cannot be directly measured in a realistic ex- 
periment. The most accessible quantity for experiments 
is the spectrum of emitted radiation. In order to eval- 
uate such spectra, one needs to account for absorption 
and phase-matching Q , which is beyond of scope of this 
work. Nevertheless, important conclusions about these 
spectra can be made by analyzing the Fourier transform 
of the polarization response P(t) that we evaluate with 
our model. In Fig. [5j we plot polarization spectra for the 
same simulations as those presented in Fig. 0J We see 
that the spectra are also sensitive to the carrier-envelope 
phase of the laser pulse. Similarly to excitation probabil- 



ities, the CEP-dependence is very weak for cjl = w g/4, 
but it becomes very pronounced already at = <^g/6, 
and it affects not only the quantum-beat part of the po- 
larization response, but also the conventional harmonics 
of third, fifth etc orders. 

Consistently with the results shown in Fig. [31 the 
quantum-beat signal rapidly decreases with decreasing 
wl- At the same time, according to Fig. 0] the excitation 
probabilities are comparable in all the four cases. Conse- 
quently, the decrease in the polarization response is due 
to dephasing. Most importantly, we used longer pulses 
for longer wavelengths in order to have the same number 
of optical cycles in all the pulses. As a result, dephas- 
ing played a more significant role for the 14-fs pulse with 
wl = Wg/10 than for the 5.5-fs pulse with wl = w g/4- In 
addition to that, the dephasing time is inversely propor- 
tional to the spectral width of the excitation, and excited 
electrons occupy a larger part of the Brillouin zone in the 
tunnelling regime (see Fig. 2]). 

IV. CONCLUSIONS AND OUTLOOK 

We have analyzed the polarization response of a model 
dielectric resembling SiC>2 to few-cycle laser pulses that 
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are strong enough to efficiently excite electrons from va- 
lence to conduction bands. These interband transitions 
create coherent superpositions of states that manifest 
themselves as quantum beats in the polarization. Ob- 
viously, the dynamics of interband excitations and de- 
phasing determine the properties of the quantum-beat 
part of the polarization response. It is less obvious that 
the instantaneous frequency of quantum-beat oscillations 
changes with time, approaching the bandgap frequency 
as the laser field attenuates. We find that, during the 
laser pulse, the quantum-beat frequency (averaged over 
a laser cycle) is larger than the bandgap frequency by 
Up/H, unless the ponderomotive energy U p exceeds the 
width of the lowest conduction band. Thus, the emitted 
burst of high-frequency radiation is negatively chirped. 
We have also observed that quantum beats appear in the 
polarization response significantly delayed with respect 
to the peak of a laser pulse, even though most electron- 
hole pairs must appear when the laser field is strongest. 
Apparently, the effect of a strong field on charge carri- 
ers suppresses single-photon transitions from conduction- 
to valence-band states. A satisfactory explanation of this 
phenomenon probably requires the development of an an- 
alytical theory, which may be triggered by this work. 

For a dielectric like SiC>2, the onset of interband tun- 



nelling (7 ~ 1) coincides with the onset of Bloch os- 
cillations (wb ~ tt^l). In the tunnelling regime, we 
have observed that the non-linear polarization and the 
final distribution of charge carriers are sensitive to the 
carrier-envelope phase of a laser pulse. In particular, 
the spectrum of the emitted high-frequency radiation 
depends on the CEP, which should be possible to ob- 
serve in experiments. Such a measurement would have 
to overcome obstacles related to the strong absorption 
above the bandgap and the phase mismatch below the 
bandgap Q. Also, it must be mentioned that the coher- 
ent few-femtosecond pulse radiated by quantum beats 
spectrally overlaps with the incoherent radiation due to 
fluorescence, which our model does not account for. Some 
ultrashort temporal gating may therefore be necessary to 
isolate the radiation due to quantum beats. Still, in spite 
of all these obstacles, such measurements should be fea- 
sible at the current stage of technology. 
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